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Abstract 

By tracking tracer particles at high speeds and for long times, we 
study the geometric statistics of Lagrangian trajectories in an intensely 
turbulent laboratory flow. In particular, we consider the distinction be- 
tween the displacement of particles from their initial positions and the 
total distance they travel. The difference of these two quantities shows 
power-law scaling in the inertial range. By comparing them with simu- 
lations of a chaotic but non-turbulent flow and a Lagrangian Stochastic 
model, we suggest that our results are a signature of turbulence. 



1 Introduction 

Understanding turbulent flow remains a challenge for both theory and ex- 
periment. With strongly coupled active degrees of freedom spanning orders 
of magnitude in space and time, intense turbulence has so far resisted a 
complete description. Although dimensional models have provided valuable 
guidance [12], we still lack a fully predictive theory of turbulence. 

Recent years, however, have seen the exciting development of new tools — 
both experimental and numerical — for measuring the Lagrangian trajec- 
tories of fluid elements in turbulence [25]. This capability has in turn 
spurred the formulation of theoretical ideas that take advantage of the 
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newly available information [llj. In particular, a broad class of geomet- 
ric models that explicitly consider the Lagrangian dynamics of the flow has 
arisen. Many of these models describe the shape of clusters of fluid elements 
[3 HI IS [28] , in an extension of classical theories of turbulent relative dis- 
persion [201 [21 [151 El US] ■ The geometric properties of individual trajectories, 
however, have also been recently been studied, viewing the trajectories as 
three-dimensional space curves and focusing on their curvature [HI [271 [8] . 

Here, we study individual particle trajectories from a different stand- 
point, blending elements of turbulent transport with geometric models. Go- 
ing back to the seminal work of Taylor [23] . single-particle turbulent trans- 
port theories have described the evolution of R(i), the vector giving the 
displacement of a particle from its position at some reference time. The 
mean-squared displacement (-R^(t)) (where R is the modulus of R) is ex- 
pected to scale as for short times and then as t for t ^ Tl, where is 
the integral time scale [19], similar to molecular diffusion processes, though 
with a flow-dependent coefficient. Displacement, however, is not the only 
distance measure we can apply to a particle path; since we follow individual 
particles in three dimensions, we can also study the arc length of the trajec- 
tories, i.e., the total distance travelled by individual particles. We denote 
this arc length by S{t). S and R are identical for simple unidirectional flow; 
in highly rotational flow, however, as is the case in turbulence, we expect 
significant differences between the two. Since coherent structures in turbu- 
lence are expected to be strongly vortical, the difference between S and R 
may possibly mark these structures. 

We find that although the statistics of S alone are similar to those of R, 
the difference of the two shows very different behavior. We measured both 
{S'^ — R'^) and {{S — R)"^) , which differ by the non-negligible covariance (SR). 
Each of these differences shows power-law behavior in the inertial range, 
though with different exponents; we find that {S'^ — R'^) ~ (similar to clas- 
sical Richardson dispersion for particle pairs) and {{S—R)'^) ~ ^3.7±o.2 
effort to isolate the contributions of turbulence to our findings (rather than 
simple kinematic effects [14j), we also measured these statistics in two syn- 
thetic flows: the Arnold-Beltrami-Childress (ABC) flow [10], which shows 
chaotic particle motion but is not turbulent, and Sawford's second-order 
Lagrangian stochastic model [22] , which approximates non-intermittent tur- 
bulent transport but does not have spatial structure. Comparing our exper- 
imental results with these synthetic flow fields, we conjecture that a large 
separation of time scales, and therefore a (Lagrangian) inertial range, is 
responsible for the detailed results we find. 

In Sec. [2] below, we briefly describe our experimental flow and the mea- 
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surement techniques we used to gather our data. We then present and 
discuss the statistics in Sec. [3l and finally summarize our findings in Sec. |H 

2 Experimental Details 

To generate fully developed turbulence in the laboratory, we used a von 
Karman swirling fiow of water between counter-rotating baffled disks. De- 
tails of the apparatus are given elsewhere [18]. By rotating the disks at 
frequencies up to 5 Hz, we drove turbulence at Taylor-scale Reynolds num- 
bers Rx = yJlbu'L/v, where u' is the root-mean-square velocity, L is the 
integral length scale, and v is the kinematic viscosity, up to 815. Here, we 
show data for R\ = 690 and 815. Cooling water circulated through the 
top and bottom plates of the apparatus in order to maintain a constant 
temperature in the working fluid. 

We measured the turbulence by tracking the motion of small, solid tracer 
particles seeded in the flow. These polystyrene particles have a density of 
1.06 g/cm^ and a diameter of 25 //m, and have previously been shown to 
be passive tracers in this flow |25] . With RMS speeds of order 1 m/s and 
accelerations of order 10^ m/s^, very fast cameras are required to record 
the trajectories accurately. We imaged the particles at rates of up to 27 000 
frames per second. We used three cameras arranged around the mid-plane of 
the tank in order to measure the three-dimensional positions of the particles 
by stereomatching the individual images. In order to study the long-time 
transport of the particles, we illuminated a measurement volume of 5 x 5 x 5 
cm^ with two pulsed, frequency-doubled Nd:YAG lasers delivering roughly 
150 W of laser light. For comparison, the integral length scale L is 7.1 cm 
in this flow. 

We used particle-tracking methods to make quantitative measurements. 
We process the raw movies with a multi-frame predictive tracking algorithm 
|17) that reconstructs the long-time trajectories of individual particles. In 
order to keep the stereomatching ambiguities at a manageable level, the seed- 
ing density of tracers is kept low (roughly 100 per image). Even so, particles 
can occasionally drop out of view, leading to gaps in the raw trajectories. 
In a post-processing step, we bridge these gaps by re-tracking the data in a 
six-dimensional space spanned by the particle positions and velocities [26j . 
In this way, we obtain longer tracks that give us access to later-time statis- 
tics. Velocities and accelerations are subsequently computed by convolving 
the measured trajectories with a smoothing and differentiating kernel [15j . 
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Figure 1: Evolution of the mean-squared displacement (i?^) and distance 
travelled {S'^) as a function of time for two different Reynolds numbers. 
Both scale ballistically (-^ t^) at short times, and trend downwards towards 
diffusive scaling (~ t) at longer times. The dashed line is a t"^ power law as 
a guide for the eye. 

3 Results and Discussion 
3.1 Distance and Displacement 

The simplest transport quantity that can be computed for a single fluid 
element is its mean-squared displacement from some initial position. This 
quantity, which we denote as (i?^), is simple to model in turbulence [23} [T9l 
|9]: we expect it to scale as n'^t'^ for t <^ Tl, where Tl is the integral time 
scale, and as u''^TLt for t ^ T^. Our measurements of {R'^) are shown in 
Fig. [H and show no unusual behavior. 

Displacement, however, is not the only way to measure the change in 
position of a fluid element. Since we know its full trajectory, we can also 
measure the total distance travelled in a time t by measuring the arc length 
of the trajectory. We label this positive-definite quantity S, and study its 
mean-square value {S"^), as shown in Fig.[T] Like (i?^), we find that (S"^) 
for short times, with nearly the same scaling constant. This behavior should 
be expected: at short times, distance and displacement are very similar. 
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Figure 2: The difference (5^ - R^) as a function of time, plotted in (a) 
logarithmic coordinates and (b) compensated by et^, where e is the energy 
dissipation rate. In the inertial range, this difference appears to scale as et^ 
with a scaling constant of roughly 2 and a scaling region that grows with 
Reynolds number. 



3.2 Difference Statistics 

More intriguing than the statistics of {Fi?) and (S^) alone are the statistics of 
their difference, for two main reasons. First, since both (i?^) and (S^) scale 
like f^ for t <^Tl, we may be able to observe sub-leading contributions to 
their scaling by subtracting them. Additionally, large discrepancies between 
the displacement of a fluid element and the distance it has travelled may be 
associated with highly vortical regions, and may be indicative of coherent 
structures such as vortex tubes. 

In Fig. 2(a) we show the simple difference (S^) — (i?^) as a function of 



time. At the two Reynolds numbers shown, this difference scales as a power 
law in the inertial range. Unlike (i?^) and (5^) alone, however, this difference 
scales as rather than t^, similar to the prediction for Richardson scaling of 
the mean-squared separation of particle pairs in turbulence. This observed 
scaling can be rationalized by applying Kolmogorov-style reasoning: {S"^) — 
{B?) has dimensions of length squared, and thus must scale like et^ (where e 
is the typical energy dissipation rate per unit mass) in the inertial range. We 
check this scaling prediction by compensating [S"^) — {R'^) by et^, as shown 



in Fig. 2(b) To check the robustness of this power law, we performed a 

(1) 



local-slope analysis by computing the logarithmic derivative 

dlog(S2 -i?2) 



dlogt 
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Figure 3: Local slopes (i.e., the logarithmic derivative) of (S*^ 
function of time, confirming the scaling in the inertial range. 



as a 



which gives the local scaling exponent of the data as a function of time. 
Robust power laws will have long plateaus in a local-slopes plot. Our results 
are shown in Fig. [3l Local slopes have been shown to be a clean way of 
identifying power laws in experimental data sets without the need to fit to 
potentially small domains [3l [1] . This analysis confirms the existence of a 
^3 j-ggi]^g_ The t"^ regime therefore appears to be a distinctive feature of 
high-Reynolds-number turbulence. 

But even though these data show a scaling regime compatible with Kol- 
mogorov's 1941 theory [12j (commonly known as K41), it is not clear that 
K41 can be applied. In K41, scaling predictions can only be made for 
scale-local quantities. Differences can often be assumed to be scale- local 
p], but this heuristic argument cannot distinguish between (5^) — (i?^) and 
((S — i?)^), two quantities that are dimensionally equivalent. The distinction 
between these two statistics is the covariance (SR); cross terms such as this 
one have been shown to play an important role in obtaining clean scaling in 
turbulent relative dispersion [Sj [18] . In Fig. [H we show our measurements 
of {SR). As expected, this covariance is not at all negligible (since S and 
R are highly correlated), and grows as a power law in time. As shown by 



a local-slope analysis in Fig. 4(b) (SR) scales as at short times and t at 
longer times, just like (5^) oi{W). 

We therefore consider the quantity {{S — R)"^), which we show in Fig. [5l 
For inertial-range times it scales as t3.7±o.2^ which is seen clearly from a 
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Figure 4: (a) The covariance (SR) plotted in logarithmic coordinates. Just 
as for (5^) and (i?^), (SR) scales as at short times and subsequently as t 
at longer times, (b) Local slopes for (SR), confirming the scaling laws. 




Figure 5: (a) The difference {{S — R)'^) plotted in logarithmic coordinates. 
The straight line is a t^'^ power law as a guide for the eye. (b) Local slopes 
for {{S — R)"^), suggesting that the inertial-range t^''' scaling is robust. 
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local-slope analysis, shown in Fig. |5(b)[ The local-slope analysis also sug- 
gests that (as expected) the temporal extent of the inertial-range power law 
scaling region grows with Rx. A scaling exponent of 3.7 cannot be obtained 
from K41; and since these statistics are only second order, any possible in- 
termittency corrections are likely to be small. As we confirm below, this t^ '^ 
regime appears to be unique to turbulence and requires the existence of an 
extended inertial range, but does not need intermittent particle accelerations 
or velocity increments. 



3.3 Comparison with Other Flows 

In order to ensure that the statistical signatures we see are truly related to 
the turbulence and not simply to kinematics [14] . we also simulated par- 
ticle trajectories in a chaotic but non-turbulent flow and in a Lagrangian 
stochastic model. 



3.3.1 ABC Flow 

The steady Arnold-Beltrami-Childress (ABC) flow is a canonical example 
of a flow field that has a simple spatial structure but can show complex, 
chaotic Lagrangian dynamics [lOJ- The velocity field, which is a solution of 
the Euler equations, is given by 

u{x, y, z) = {A sin kz + C cos ky) 

+ {B sin kx + A cos kz) ey 

+ {C sin ky + B cos kx)ez, (2) 

where A, B, and C are constants, k is the spatial wavenumber of the flow, 
and e^;, e^, and are unit vectors in the three Cartesian directions. Even 
though the spatially periodic flow field is steady, fluid elements advected in 
this flow will decorrelate on a time scale related to the time necessary to 
exit a unit cell. We estimate this time scale to be 

2Tr/k 

""^'^ = [A^ + B^ + C^f-- 

We simulated a population of fluid elements evolving in the ABC flow 
using a fourth-order Runge-Kutta algorithm, choosing A = 12, B = 15, and 
C = 10 [2Tj, and measured S{t) and R{t) along their trajectories. In Fig. [6l 
we show the local slopes for {(S — R)'^) and (5^ — i?^) for these fluid elements. 
At long times, both quantities saturate to a scaling. This result is to be 
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expected since for t » Tabc, (R^it)) ~ t and {S^{t)) ~ t^; at long times, 
the leading-order t'^ term dominates. For shorter times, we observe no clean 
power-law scaling. Each curve in Fig.[6]has a short region that may indicate 
a short power-law regime; the exponents, however, are very different from 
the and t'^'^ scalings seen in the turbulence data. The ABC flow misses 
some of the ingredients of turbulence. 

3.3.2 Lagrangian Stochastic Model 

Although our measurements in the ABC flow do not show the same scaling 
as the turbulence measurements for (S^ — i?^) and {{S — R)^), it remains to 
be seen whether our observations are a result of the broad spectral activity 
in turbulence (that is, the presence of both a large and a small time scale 
separated by an inertial range) or are due to intermittency. We therefore 
used Sawford's second-order Lagrangian stochastic model [22] to test these 
possibilities. 

Sawford's model assumes that the material derivative of the acceleration 
is a white-noise process, so that the Lagrangian acceleration is given by a 
stochastic differential equation. The velocity and position of a fluid element 
are obtained by successive integrations of the acceleration. Though this 
model has two time scales (corresponding to the integral time and the 
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Figure 7: Local slopes for {S'^—R'^) and {{S—R)"^) computed using Sawford's 
model [22j . The observed turbulent inertial-range scaling is evident here, 
despite the absence of intermittency in the model. 

Kolmogorov time r^) and therefore has an inertial range, it is not intermit- 
tent: the probability density functions of the acceleration and the velocity 
increments are Gaussian. 

Our results computed using this model with Rx = 815 (averaged over 
many realizations) are shown in Fig. [71 We see behavior much more similar 
to the turbulence than to the ABC flow, though with shorter scaling ranges 
at similar Reynolds numbers. These results suggest that we are observing 
normal inertial range scaling rather than the anomalous scaling associated 
with intermittency, a result that is to be expected for low-order statistical 
moments like the ones we consider here. 

4 Summary 

We measured statistics of the trajectories of single particles in an intensely 
turbulent laboratory flow, focusing on the displacement of the particles from 
their initial positions and the total distance travelled by the particles. The 
mean-squared values of each of these lengths scales as expected, growing 
ballistically at short times and diffusively at long times. Differences of the 
distance and displacement, however, reveal new scaling laws that depend 
on exactly how the difference is taken. Although we do not have a detailed 
theoretical explanation for our observations at this time, we conjecture that 
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our results are a signature of the separation of time scales present in tur- 
bulent flow, and may possibly be signatures of Lagrangian structures in the 
flow. We hope that these results stimulate theoretical investigations of the 
measured statistics. 
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